Two-band BCS model describes well the thermodynamics of MgB 2 
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Based on weak-coupling anisotropic BCS theory, the temperature dependence of energy gap and 
the specific heat are evaluated for MgB2 superconductor, and the results are compared with exper- 
imental data. We show that the weak-coupling anisotropic BCS theory describes thermodynamic 
experimental data with high precision 3-6%. 
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A keen interest excited by discovery and experimental 
investigation of a new high T c superconductor MgB 2 is to 
a large extent associated with its dissimilarity to cuprate 
superconductors. The superconductivity of MgB 2 is def- 
initely three-dimensional effect, whereas in cuprates it is 
presumably two-dimensional (2D). Nevertheless, the su- 
perconducting gap in MgB2 displays strong anisotropy. 
The most precise tunneling measurements by Gonnelli et 
ai. give the value 2.6 for the ratio of the gaps at two 
conductivity bands. On the other hand, the measured 
gaps are the same for the tunneling in afc-plane and in 
c-direction, indicating that they do not depend on direc- 
tion within each piece of the Fermi surface. 

An important problem is how strong is the interac- 
tion in MgB 2 . First-principles calculations Q-EI DI indi- 
cate that electron-phonon interaction is not weak and 
that Eliashberg description is appropriate. However, 
anisotropy and interaction were shown to influence ther- 
modynamics oppositely. For example, the anisotropy de- 
creases the relative discontinuity of the specific heat at 
the transition point 0, 0], whereas the first correction 
due to interaction increases it Q. Besides, MgB2 is a 
very hard material with high value of Debye frequency, 
which usually correlates with a weak coupling. There- 
fore it is not a priori clear what is more substantial in 
the case of MgB 2 . 

The purpose of our work is to demonstrate that the 
anisotropy effects are more substantial at least for ther- 
modynamic measurements. We show that, as a matter 
of fact, the weak coupling anisotropic BCS theory de- 
scribes all known thermodynamic experimental data in- 
cluding the temperature dependence of the energy gap 
and specific heat with a high precision 3-6%. 

The main features of anisotropic weak coupling BCS 
model were elucidated in the early 1960s |a, la, S [' 
the ultimate result being the factorization of the gap 
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A(T,k)=Q(T) X (k), (1) 

which was experimentally verified by Zavaritskii 
The function of angle x(k) is the eigenfunction of the 



interaction operator V(k, k') corresponding to the max- 
imal eigenvalue A + . ft satisfies linear homogeneous inte- 
gral equation: 



V(k,k') X (k') 



da' 

VpVF 



A+X(k). 
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Integration in Eq. proceeds over the Fermi surface 
with da — and dS being a differential area of the 
Fermi surface; Vp = J is the electron density of the 
state per spin at Fermi level. The function x(k) is nor- 
malized as follows: 
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The angular average value (X) is: (X) = J -^T"- the 
temperature dependent factor Q(T) can be found from 
the orthogonality condition: 
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The value Q(0) is associated with the transition temper- 
ature T c by the following relationship: 

^ = -cxp(-( X 2 (k)ln| X (k)|)), (6) 

here 7 — e c — 1.781072 • • • and C is Euler's constants. 
The specific heat C(T) reads: 



(7) 



where G(x) = 2x J °° cosh(2<p) F(x cosh ip) dtp. 

We now apply these formulas to MgB 2 . The Fermi 
surface of MgB 2 has two cr-type 2D cylindrical hole 
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sheets and two 7r-type three-dimensional tubular net- 
works 0, 0, 0] . We accept a simple model introduced 
first by Moskalenko 0] , in which the interaction does not 
depend on the momentum inside each band, but only on 
the band index. Thus, it can be written as 2 x 2 Hcrmitian 
matrix Vik (i,k — a, n). The order parameter (energy 
gap) in each band in such a model does not depend either 
on the momentum within each band and can be described 
by a 2D vector with components A^, A w . The validity 
of this simple model is supported by the tunneling mea- 
surements of the energy gap 1] . which displays the same 
values for two gaps in afo-plane and in c-direction. The 
normalized wave function of the Cooper pairs Xk nas the 
same property: Xo-(k) = X<r, 36r(k) = Xn, where x<r and 
Xtt are two constants. We introduce an additional sim- 
plification assuming these constants to be real. Let us 
denote the density of states in the a and tt bands as vpa 
and Vf-ki respectively. Then the definition of an aver- 
age value (X) for any physical value X, which does not 
change within each band reads: 



(X) = X a c a + 



(8) 



where c a and c w are statistical weights of the bands c a — 
vfoIvf and = vfw/vf, vp = vp a + vp^. The general 
normalization condition Eq. Jjyl for this model reads: 
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+ X-K C * = 1 ■ 

Equation © can be written explicitly as follows: 

0(0) 7T 



where % av = X%° C *X%* C 
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We assume the values c a 
().-',-", and = 0.56 as found from density-functional the- 
ory calculations in Refs. [121 Il5| . The second fitting 
parameter is T c . There is no experimental discrepancy on 
this value, and it is commonly accepted to be T c ss 39 K. 
One additional fitting parameter for the two-band the- 
ory is the ratio S = Xo/x-n- We have extracted it from 
the tunneling gap measurements extrapolating them 
to zero temperature: 



5 = Xa/Xn ~ 2.54. 



(11) 



Equations l|10|l and allow us to determine x<? and 
X-k separately : x<j = S/y/c a 5 2 + = 1.38; x-n = 
1/^/caS 2 + c> ~ 0.54. According to the weak-coupling 
theory, the ratio 5 must be the same at any tempera- 
ture. This crucial condition is satisfied in the tunneling 
experiment Q with all experimental precision. 

For the temperature dependence of the gap in the BCS 
two-band model, we find from Eq. (0J: 
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FIG. 1: The solid curve depicts the ratio Q(T)/Q(0) vs t = 
T/T c for the two-band model; the dashed curve is the same 
value for the standard (isotropic) BCS theory. 
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FIG. 2: The solid curve is the theoretical graph of A CT vs 
T /T c ; the dashed curve is the same for A^; "+" and "x" 
represent experimental data by Gonnelli et al. Q 



Here q{t) = Q(t)/Q(0) and t = T/T c . F(x) is defined by 
Eq. ©. The graph of the function q(t) is shown in Fig.Q] 
by the solid curve. The dashed curve in Fig.Q]represents 
q(t) in the isotropic single-gap model (standard BCS 
model). The graphs of the energy gaps A a = Q{t)xa 
and A,,. = Q(t)xrr vs T/T c are shown in Fig. [21 together 
with the experimental data 0, which agree with theory 
within the limits of experimental uncertainty. 

The specific heat in the two-band model is given by 
the following equation directly stemming from Eq. 0: 



C(T) 
C N (T) 



12 [cvxlr a (y„) + c n xlr a (y^)} 2 

7C(3) c a xin(ya) + CnX%r b (y w ) 



(13) 
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FIG. 3: The solid curve is the theoretical graph of the specific 
heat for the two band MgEb vs t — T/T c ; the circles are the 
experimental data due to Bouquet et al. the dashed curve 
is the theoretical plot of the specific heat given by the isotropic 
BCS theory. 



where Cn(T) = 7T is the specific heat for the normal 
metal; ^ = f^f^>^ = is : f^- Tlie functions n are 
defined by integrals r t (x) = f_™ gi{\J x 1 + y 2 )dy : i = a, 
b, c, where gi read: 
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7T 2 COsh 2 ; 



(14) 



For technical details related to this calculation see Mis- 
honov et al. |l7j : the functio ns q ,; were introduced and 
graphically presented in R ef. 1 181 The jump of the spe- 
cific heat at T c reads: cf.|5lll4| 



AC(T C ) 12 (x 2 c* + xlc*) 



C N {T C ) 7C(3) 



(15) 



For the data specified earlier, we find AC(T c )/Cn(T c ) = 
0.874. It agrees with the high precision measurements by 
Bouquet et al. 0] with about 3% precision. In Fig.|3]the 
ratio C(T)/C N (T) vs. T/T c is plotted. The solid curve 
is the prediction of the two-band weak coupling theory; 
the dots are experimental data by Bouquet et al. [l9| . 
courteously sent to us by the authors. The theoretical 
graph C(T)/Cn vs T/T c agrees well with the experimen- 
tal data everywhere except of a range of low temperature 
T/T c < 0.2. The discrepancy most probably is caused by 
a relatively small variation of the gap within one band. 
The specific heat at low temperature is proportional to 
e -A mi „/T^ wnereas tunneling measurements give the 
value of the gap along the direction of the tunneling. 



Given the value of discrepancy, we can estimate the vari- 
ation of the gap A- A min ~ 0.1-0.15 T c In 2 « 3.3-4.2 K. 
It is about 8-12% of the value of the smaller gap. 

Another group of available experimental thermody- 
namic data relates to magnetic properties: the energy 
gaps in external magnetic field |20j and the dependence 
of the second critical field on temperature [2l| . The de- 
pendence of H C 2 on temperature was considered theoret- 
ically in the framework of anisotropic BCS model by two 
groups of authors |22j, |23j . based on classical approach 
by Helfand and Wertheimer Unfortunately, a con- 

sistent solution of these problems at any temperature be- 
tween and T c requires much more detailed knowledge 
about the Fermi surface. For example, to reach a sat- 
isfactory convergence Miranovic et al. |22l | were forced 
to introduce 11 different parameters characterizing the 
Fermi surface and electron interaction. It is clear, that 
our real knowledge of the Fermi surface is too poor for 
such a sophistication. Dahm and Schopohl |23| applied 
a simplified model of the Fermi surface as consisting of 
a torus and cylinder characterized by 4 parameters only 
and assumed a plausible variational procedure introduc- 
ing one more parameter. As it could be expected from 
the results by Miranovic et al., the number of parameters 
is too small to ensure a reasonable precision. Indeed, a 
satisfactory agreement with the experiment in Ref. |23j 
is reached at the expense of a rather exotic choice of pa- 
rameter. Summing up, the magnetic properties can not 
be described by such an elementary theory as the de- 
scribed above two-band BCS model and require much 
more sophisticated approach even in the weak coupling 
approximation. 

Let us discuss why this simplified theory works so well. 
Let us start from the assumption supported by experi- 
ments that the gap does not vary within each band. The 
in-band isotropy of the gap could be a result of suffi- 
ciently strong in-band scattering. At the scattering time 
t ~ 10~ 14 s, i.e. at the residual resistance larger than 
10 -5 ficm, the energy gap becomes isotropic. However, 
the ratio of the gaps for different bands still remains big- 
ger than 2 indicating that the inter-band scattering must 
be much weaker. It should be emphasized that it is the 
density of states which becomes isotropic, whereas the or- 
der parameter remains anisotropic unless the Ioffe-Regel 
limit tef ~ 1 of scattering rate is reached |25j. The 
tunnelling experiment measures just the density of state. 

Second question is why the weak-coupling model gives 
so high accuracy. Two different aspects must be enlight- 
ened. First, the separability of variables for the order 
parameter, even in the framework of the weak-coupling 
approximation, has the precision of of the weak coupling 
constant, i.e. (ln-^-) _1 ~ 0.3. For the case of the two- 
band model such a crude estimate can be checked more 
accurately by a direct solution of the nonlinear matrix 
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equation for the energy gap. It has a following form: 
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where i,j 

tanh u 



take 



tanh \/ n 2 +z 2 



values cr, 7r and /(a;) = 
c?u. Its solution can be found 



as a superposition of two normalized eigenstates of the 
corresponding linear equation: Aj = Q+VP+j + Q_^_j. 
In our calculations we used only one of them, ty+ 
corresponding to the larger eigenvalue A+. Such an 
approximation is justified when the second eigenvalue 
A_ is much less than A+, even if A+ is not very 
small. Indeed, the symmetrized matrix V with ma- 
trix elements Vij — ^CiCjVij can be represented as 
V = A+|+)(+| + A_|— )(— |. This representation shows 
that, at A_ = the operator V is separable, and the 
solution of non-linear equation i|16[l is factorizable: 
A = Q(T)^/ + . The equation A_ = is equivalent 
to DetV = VcrcrV^ - V£t = 0. Though such a fine 
tuning of parameters seems improbable, our numerical 
calculations demonstrate that the ratio A_/A+ and the 
thermal variation of of the ratio A W /A CT remain small 
(about 3%) even at DetV / (c„V a(7 + c^V^) ~ ±0.2. 
Thus, the experimental facts seem to indicate that one of 
the two eigenvalues is significantly smaller than another. 
Such a situation occurred earlier in a band calculation 
for high-T c -superconductors |27j . 

The second aspect mentioned in the preamble is that 
the BCS approximation itself has a low precision and 
should be substituted by the Eliashberg formalism. The 
numerical calculations by Golubov et al. Q indicate that 
the Eliashberg weight functions are very small in a broad 
range of low energy and has rather sharp peaks in the 
range of 800-1000 K. This is an unusual situation. Leav- 
ens and Carbotte [28| considered an extended Eliashberg 
weight function a 2 F(u) centered at values u ~ ujq much 
larger than the superconducting energy gap A(0). They 
argued on the basis of numerical calculations that in this 
case the function A(u>) varies very weakly at u < luq and 
then rapidly changes sign. They even modeled A(w) by 
the step function. Their arguments seem to be correct 
for the considered case as well. Then it is obvious that by 
integrating in the range of high frequency, it is possible to 
obtain the BCS-like equations with a renormalized, not 
small interaction between electrons with momenta on the 
Fermi surface. Though such an explanation is plausible, 
further study of the Eliashberg equation with a model 
weight is highly desirable. 
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